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ABSTRACT 

The  existence  of  a non-negative  solution  of  the  boundary-value 
problem 

y"  + fy'  + [y  - d+4)]y  = o , 

r 

y(0)  = 0,  y(oo)  = o , 

is  proved  by  a shooting  argument.  The  equation  arises  in  an 
electromagnetic  theory  for  strong  interaction  in  charged  media. 
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THE  EXISTENCE  OF  A NON-NEGATIVE  SOLUTION  OF  AN  ORDINARY 
DIFFERENTIAL  EQUATION  ARISING  IN  ELECTROMAGNETIC  THEORY 


J.  B.  Me  Leod 

1.  Introduction 

We  are  interested  in  the  boundary -value  problem  associated  with 

y"  + j y‘  + [y  - (i  + \ )]y  = o (*•  U 

r 

where  the  function  y(r)  satisfies  the  boundary  conditions 

y(0)  = 0,  y(»)  = 0 . (1.  2) 

This  problem  arises  in  an  electromagnetic  theory  of  strong  interaction  which 

is  discussed  by  Bergstrom  in  [1],  although  the  equation  (1.1)  is  not  specifically 

derived  in  [1],  and  I am  indebted  to  Seymour  Parter  for  bringing  the  problem  to 

my  attention.  We  want  to  show  the  existence  of  a non-negative  solution  to 

this  problem,  and  the  method  of  attack  is  to  "shoot"  from  the  origin.  To  this 

end  we  prove  the  following  sequence  of  lemmas,  which  we  state  now  with  the 

proofs  given  later. 


! 


Lemma  1.  Given  any  a > 0,  there  exists  a unique  solution  of  (1.  1)  which 
has  the  asymptotic  behaviour 

y(r)  ~ ar  as  r — 0 , 

and  plainly  this  solution  is  non-negative  initially  (i.  e.  for  r sufficiently 
small). 

Lemma  2.  Any  solution  y(r)  o_f  (1.  1)  which  is  initially  non-negative  must 
possess  one  of  three  mutually  exclusive  properties: 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024. 
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(i)  it  remains  non-negative  for  only  a finite  range  of  values  of  r ; 


I 

t 


(ii)  it  remains  non -negative  for  all  r,  and  y(r)  — 1 as  r 

(iii)  it  remains  non-negative  for  all  r,  and  y(r)—  0 as  r 


00 


— m?la  l • The  set  of  a(  >0)  for  which  the  solution  given  by  Lemma  1 has 
the  property  (i)  of  Lemma  2 is  denoted  by  S,  and  is  open  in  the  topology  of 
the  positive  semi-axis. 

— mma  4-  The  set  of  a (>0)  for  which  the  solution  given  by  Lemma  1 has 
the  property  (ii)  of  Lemma  2 is  denoted  by  and  is  open  in  the  topology 
of  the  positive  semi-axis. 


^g-mma  5v  II  Q'(>0)  is  sufficiently  large,  then  the  solution  given  by 
Lemma  1 has  the  property  (i)  of  Lemma  2,  so  that  e S . 

Lgmma  6-  II  «(>0)  is  sufficiently  small,  then  the  solution  given  by 
Lemma  1 has  the  property  (ii)  of  Lemma  2,  so  that  a € S2  . 


Once  Lemmas  1-6  are  proved,  as  they  are  below,  then  the  following 


existence  theorem  becomes  an  easy  consequence. 


H?g°reiin-  There  exists  a non-negative  solution  to  the  boundary -value  problem 
(1.  1 ) — ( 1 - 2),  the  solution  being  in  fact  strictly  positive  except  at  r = 0 . 

^ro°^'  Lemmas  3-6  imply  that  the  sets  S^,  are  open  and  non-empty. 

Since  the  positive  semi-axis  is  connected,  not  all  a > 0 can  belong  to 
Sj,  S2  . Hence  there  exists  at  least  one  value  of  a which  is  in  neither 

Si  nor  S and  for  which  therefore  the  corresponding  solution  given  by 
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Lemma  1 has  the  property  (iii)  of  Lemma  2.  This  completes  the  proof,  ex- 
cept for  the  trivial  observation  that  the  solution  must  be  strictly  positive 


- 

Sf 

.* 


•f 


for  r > 0,  since  if  it  were  not  it  would  have  to  touch  the  value  0,  at 
r = rQ,  say,  when  y(rQ)  = y'(rg)  = and  this  implies  from  (1.1)  that 
y = 0,  contradicting  a > 0 . 

It  is  perhaps  worth  remarking  that  there  is  a difference  in  character 
between  Lemmas  1-4  on  the  one  hand  and  Lemmas  5-6  on  the  other.  Thus 
Lemma  1 is  a statement  about  the  asymptotic  behaviour  of  solutions  of 
the  equation  at  the  origin,  and  this  is  something  that  an  analyst  can  ex- 
pect to  discuss  relatively  easily,  even  for  a much  more  complicated 
equation  than  (1. 1),  provided  only  that  the  coefficients  in  the  equation  are, 
asymptotically  at  the  origin,  of  a fairly  simple  character.  In  the  same 
way,  Lemma  2 is  concerned  with  the  asymptotic  behaviour  at  infinity, 
and  Lemmas  3-4  are  statements  about  the  continuity  of  the  solution  in  the 
initial  parameter  a . All  these  lemmas  are  therefore  accessible  to 
analytical  proof,  again  even  for  a more  complicated  equation;  and  con- 
versely, since  Lemma  2 is  a statement  about  all  (and  so  an  infinite 
number  of)  solutions  of  a certain  kind,  and  since  infinity  and  continuity 
are  not  concepts  which  the  computer  recognizes,  it  is  difficult  to  see  how 
these  lemmas  could  be  established  other  than  analytically. 

The  purpose  of  Lemmas  5-6,  however,  is  only  to  show  that  there 
exists  one  solution  in  the  set  Sj  and  another  solution  in  the  set  . 
(More  than  this  is  stated  (and  proved)  in  these  lemmas,  but  the  existence 
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of  one  solution  in  each  set  is  all  that  is  required  in  the  application. ) 

To  establish  such  results,  a solution  has  to  be  followed,  not  only  for 
small  and  large  r,  but  also  in  general  for  moderate  r as  well,  and  in 
the  case  of  a complicated  equation,  this  may  well  prove  impossible  for 
the  analyst.  Even  in  the  present  case  the  analysis  is  delicate.  On  the 
other  hand,  since  we  are  concerned  with  only  a finite  number  of  solutions, 
it  may  well  be  possible  for  the  computer  to  follow  these  solutions  and, 
with  the  aid  of  interval  analysis  or  some  comparable  argument,  prove 
that  the  solutions  so  followed  satisfy  the  inequalities  that  place  them 

* 

in  the  required  sets. 

The  point  of  these  remarks  is  that  a shooting  argument  thus  breaks 
naturally  into  two  parts,  one  of  which  should  be  accessible  to  the  analyst 
and  the  other  to  the  computer.  It  becomes  therefore  a very  powerful  tool 
for  the  proof  of  existence  of  solutions  of  two-point  boundary-value  problems. 

2.  Proof  of  Lemma  1 

We  can  write  (1.  1)  in  the  form 

y"  + ~ y*  - \ v = (i  - v)y  , (2. i) 

r 

and  since  two  linearly  independent  solutions  of 

2,2 

y"  + - y - — y = 0 
r 

2 

are  y = r,  y = 1/r  , the  variation  of  constants  formula  shows  that  (2. 1) 

is  equivalent  to  the  integral  equation 

2 l r t3 

y = Ar  + B/r  + ^ f (r  - -j  ) {1  - y(t) } y(t)dt 
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where,  for  any  fixed  r^  > 0,  all  solutions  of  (2.  1)  which  exist  in  some 
neighbourhood  of  r = rQ  are  given  by  arbitrary  choices  of  the  constants 
A,  B . In  fact,  since  we  are  interested  in  solutions  for  which  y ~ ar  , 
we  can  actually  choose  r^  = 0,  since  the  integral  continues  to  have 
meaning  then,  and  we  shall  then  have  to  choose  B = 0,  A = a . Thus 
such  a solution  satisfies  the  equation 

1 / t3 

y -■  ar  + — J (r-  — ) {1  - y(t)}  y(t)dt  , (2.2) 

0 r 

and  the  usual  approach  of  solution  by  iteration  shows  that  this  equation 
has  one  and  only  one  solution  for  which  y ~ »r,  and  so  the  lemma  is  proved. 

3.  Proof  of  Lemma  2 

We  first  note  from  (2.2)  that 

1 / t3 

y < ar  + - J (r  - — ) y(t)dt  , 

0 r 

from  which  it  follows  that  y < Y,  where  Y is  the  solution  of  the  integral 


equation 


= «r  + T / (f  - ~ ) Y(t)dt  , 


with  Y ~ ar  as  r -*  0 . But  it  is  easy  to  deduce  from  (3.  1)  that  Y 
satisfies  a linear  differential  equation,  and  so  exists  for  all  r and  is 
bounded  in  any  compact  set.  Since  a non-negative  solution  y of  (2.  2) 
is  bounded  below  by  0 and  above  by  Y,  it  follows  that  y exists  for 
all  r and  is  bounded  in  any  compact  set.  Hence,  if  alternative  (i)  in 
the  statement  of  the  lemma  does  not  hold,  then  certainly  y exists  for 
all  r and  we  have  only  to  prove  that  y-*l  or  y -*  0 as  r -*  x . 


We  therefore  assume  for  the  remainder  of  the  proof  that  alternative 
(i)  does  not  hold,  and  show  first  that,  if  y is  ultimately  monotonic, 
then  y -*  1 or  y -*  0 . For  if  we  suppose  for  contradiction  that 
y(oo)  = k ^ 0,  1,  where  y(oo)  exists  since  y is  monotonic,  then  either 
0 < k < 1 or  k > 1 (with  possibly  k = +»)  . If  0 < k < 1,  (2.  1)  gives 

(r2y')'  = -r2y{y-(l  + 2/rZ)}  , 

so  that 

(r2y')'  ~ r2  k(l  - k)  , 

2 13 

r y'  ~ j r k(l  - k)  , 

y'  ~ \ rk(i  - k)  , 

which  contradicts  y -►  k,  and  a similar  argument  disposes  of  the  case 
k > 1 . 


Thus  if  y is  ultimately  monotonic,  then  y -*  1 or  y -*  o . It  re- 
mains to  discuss  the  case  where  y is  not  ultimately  monotonic.  Then 
multiplying  (2. 1)  by  y'  and  integrating  from  some  value  rQ  > 0,  we  have 


i-l  ,2,r  r y 1-1  r1  Jir  r*  £,r 

t + f -rLidt  + (iy  ir  lr  - 


r (t)  . rl  3,r  rl  2,r 

T J 

0 r_  ‘0 


. 2 j m m dt  = o . 


(3.2) 


Now 

r 

y(r)  = f y'(t)dt 
0 


by  Cauchy-Schwarz. 
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Hence 


2 / dt  = [\  / + 2 / if  dt 


ro  1 


2 1 J 3 

1 ro  ro  1 


r t 


£ [~]  +2  / (/  y'Z(u)du)^ 


ro  0 


/ „rl  .2. 


r 2 


= t^2  1 - 2ff/  y'  (UP  +2/  P 

1.  *■  A 1 /\  L 


0 r„ 


so  that  (3.2)  gives 


r1  .2 

1 2 7 + 


t y3  - t y2  - + - / y,2(u)du]  < o . 

‘ 


(3.  3) 


3 2 

Since  y dominates  y for  large  y,  this  last  inequality  plainly  implies 
bounds  on  both  y and  y‘  as  r -*  °°  . 

Now  we  are  assuming  that  y oscillates,  and  we  may  further  suppose 
that  it  oscillates  finitely,  i.  e.  that  y(°c)  does  not  exist.  For  if  y(°°)  exists, 
then  the  argument  used  when  y is  monotonic  will  show  that  either  y(°°)  - 0 
or  y(°o)  = 1 . From  the  now-established  boundedness  of  y and  y',  we  see 
from  (3.2)  that 

f dt  (3.4) 


is  bounded  as  r -»  and  since  it  is  monotonic,  it  must  converge.  We 
thus  have 


1 2 13  12 

2 y fiy  -zy 


L,  say,  as 


(3.  5) 


r 


and  three  cases  now  arise,  depending  on  whether  L = 0;  L = -1/6;  L * 0, 
-1/6  . (The  cases  L = 0,  L = -1/6  are  the  cases  in  which  the  roots  of 

13  12 

T y -2y  -L  = 0 

are  not  distinct.  ) 

If  L = -1/6,  then  the  possible  extreme  values  of  y (in  the  limit  as 
r -*00),  i.  e.  the  roots  of  the  equation 


13  12  1 

3 y - jy  *z-o  , 


are  readily  verified  to  be  1 (twice)  and  -1/2  . Since  y is  non-negative, 
this  implies  that  y - 1 . (Alternatively,  if  L = -1/6,  then  (3.  5)  gives 


^ y'2  + ^(2y  + l)(y  - 1)2  - o , 


from  which  it  immediately  follows  that  y1  -*•  0,  y -*•  1 • ) 
If  L * 0,  - 1/6,  then  the  roots  of 


1 3 1 2 

— y - — y - L = 0 

3 2 


are  distinct,  and  none  of  them  are  0 or  1 . In  the  limit  as  r-**,  y 
oscillates  between  two  of  these  roots  as  extreme  values,  and  in  the  course 
of  a complete  oscillation,  y'  is  small  if  and  only  if  y is  near  these  ex- 
treme values.  Further,  from  (1.  1),  y"  is  not  small  when  y is  near  the 
extreme  values,  since  the  extreme  values  are  neither  0 nor  1 . Hence 
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in  the  course  of  a complete  oscillation,  y'(r)  is  small  if  and  only  if  r 
is  near  a value  for  which  y'(r)  = 0,  and  so  y'(r)  is  not  small  for  most  of 
the  r-range  of  the  oscillation.  From  this  it  follows  easily  that 

ro 


is  not  bounded  as  r-*-  °o,  and  this  is  a contradiction  to  what  we  proved  at 

(3.4). 


If  L = 0,  the  possible  extreme  values  are  0 (twice)  and  3/2  . 

There  is  no  immediate  reason  (as  there  was  with  L = -1/6)  why  the  solution 
should  not  oscillate  (in  the  limits  as  r -*><»)  between  0 and  3/2,  and  the 
argument  employed  when  L ± 0,  - 1/6  does  not  apply  because  it  is  prima 

facie  possible  for  the  solution  to  remain  small  (y,  y',  y"  all  small)  for 
long  intervals  of  r,  increasingly  long  every  time  the  oscillation  is  repeated, 
and  this  situation  does  not  necessarily  contradict  the  boundedness  of  (3.4). 

To  argue  that  the  solution  cannot  oscillate  indefinitely  between  0 and 
3/2,  we  must  therefore  look  more  deeply.  Consider  the  function 


4 1 ,2  13  1 2 

F(r)  = -y*  + -y  - - y 


(3.6) 


We  are  supposing  (for  contradiction)  that  F(r)  -*  0 as  r ~»oo,  while  y 
oscillates  indefinitely  between  0 and  3/2  . As  r increases  from  a local 
maximum  rQ  of  y to  a local  minimum  r^.  so  that  y'  < 0 in  (r0>rj)>  it 
is  clear  from  (3.  2)  that 


F(rL)  - r(rQ)  < 0 . 
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A 


so  that  F decreases  between  the  maximum  and  the  minimum.  Further,  at 


a local  minimum  of  y (at  least  once  r is  sufficiently  large  that  y is  small 

12  1 3 

at  the  local  minimum),  F < 0,  since  y'  = 0 and  y dominates  — y . 

2 3 

If  therefore  we  can  show  also  that  F decreases  as  r increases  from  a 


local  minimum  of  y to  a local  maximum,  then  we  are  done,  since  F 
(evaluated  at  local  maxima  and  minima  of  y)  will  be  negative  decreasing, 
and  so  cannot  tend  to  0 . 

The  decrease  of  F between  a local  minimum  and  a local  maximum  of 
y will  be  established  from  (3.  2)  if  we  can  show  that 


f‘  . /'«Sdt  >0  , 


(3.7) 


rQ  being  the  minimum  and  r^  the  maximum.  We  start  by  demonstrating  that, 
given  <r  > 0,  small  but  fixed,  we  have 


y(r)  < Ky(s),  rQ  i s < r 1 r1  . r - s 5 <r  , (3.8) 


where  K is  a (positive-valued)  function  of  <r  but  not  of  r , r^,  i.e.  K 
is  independent  of  how  far  out  the  oscillation  is.  For  from  (1.  1)  we  have,  for 

t£  fro’ri]  ’ 

y"  (t)  < 2y(t)  , 

at  least  if  rQ  > n/T*  , which  we  may  suppose,  and  then  by  integration 


1 2 1 2 *■ 

ijr  I <fy  ) 

0 0 
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and 


y'(t)<y(tW2  . (3.9) 

Then  y(r)  - y(s)  = (r-s)  y'(£),  s < £ < r , 

< cr\/2  y(r)  , 


from  which  (3.  8)  follows. 

We  next  establish  that 

y'(t)  > Ky(t)  for  t e [rQ  + (T,  Tj  - 0-]  , (3.10) 

K again  a (positive-valued)  function  of  (r  but  not  of  rQ,  r^  . (K  may  not 
be  the  same  function  as  before,  but  it  will  cause  no  confusion  to  use  the 
same  notation  for  both.  ) For,  since  F(r)  -*  0 as  r we  can  readily 

deduce  that  y'  is  bounded  away  from  0 (and  so  a fortiori  y‘  > Ky)  so  long 
as  y is  bounded  away  from  0,  3/2  . Since  y"  is  bounded  away  from  0 
when  y is  near  3/2,  the  statement  that  y is  bounded  away  from  3/2  is 
eauivalent  to  the  statement  that  r is  bounded  away  from  r , and  so  we 
certainly  have 


y'(t)  > K y(t) 


if  t « [rQ  + o'.  rj  * o-]  and  y(t)  is  bounded  away  from  0 . But  if  tc  [rQ,  r^] 
and  y(r)  is  near  0,  say  y(T)  < , then,  from  (1.  1), 


y"(T)  > jY(t)  > ~ y(r)  - t 

^ 0 


so  that  , integrating  over  [r^,  r]  with  r > r^  + <r,  we  have 
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' 


— ■ 


1 r 2 

y'(r)  > - J y(t)dt  - — {y(r)  - v(r  ) } 
ro  ° 

1 rr  , 2 

> 2 i y(t>dt  - ~ y(r) 


1 2 
>70-  y(r-o-)  - — y(r) 

ro 


> (jx  - — ) y(r),  from  (3.  8)  , 


which  proves  the  required  result  since  cr  (and  K)  are  independent  of  r 
and  rQ  can  be  chosen  so  large  that 


2 cr 
rQ  " 2K  • 


If  we  now  consider  the  difference  in  (3.7),  we  see  that  it  exceeds 


rr<r  2 


r X yyj.  r yyl  r yyl  ^ 
J ( t ' 2 *dt  " J 2 dt  ' J 2 dt 


ro  1 


Tj  -cr  t 


rrir  . V" , ri 

>(V1}1  T~dt  • / ^d*-/ 

V"  ‘ ' ro  1 rr<r  1 


from  (3.10)  , 


ro+2(r  rr<r 


> { K( rQ  + cr ) - 1 } {/  + f } Zf  dt 


r -hr  r^-2(r  t 


ro+<r  , ri 

-/  ^d<-/ 

r_  t r, -0-  t 


(3 


Considering  sufficiently  the  integrals  involving  rQ,  we  see  that 

V2<r  , r0+2<r  2 

/ dt,  from  (3.10), 


ro+<r  1 


ro+<r  1 


> Ko- 


y (rQ+<r) 


(r0+2<r) 


2 ’ 


while 


V*  2 


J ^-dt  < ^ 2 / y-  dt,  from  (3.9)  , 


ro  1 


ro  1 


_ y <rn+(r) 

<^2  1 


from  which  it  is  clear  (for  rQ  sufficiently  large)  that  the  contribution  from 
these  integrals  to  (3.11)  is  positive.  With  similar  arguments  for  the  integrals 
involving  r^  we  see  that  (3.17),  and  so  the  lemma,  are  proved. 


4.  Proof  of  Lemma  3 

This  is  an  immediate  consequence  of  the  fact  that  the  solution  y(r,  o) 
of  (2.  2)  is  continuous  in  a . Thus  if  y(r  , Q0)  < 0,  we  must  also  have 
y(rQ,a)  <0  for  a sufficiently  close  to  . 


5.  Proof  of  Lemma  4 


Let  us  suppose  that  for  a - the  corresponding  solution  y(r,a^) 
given  by  Lemma  1 has  the  property  that  y(r,  o^)  -*  1 as  r -*  <*  . We  have  to 
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show  that,  for  a su 


r 


show  that,  for  a su  ently  close  to  o^,  we  still  have  y(r,„)  -*•  1 as 

r — or  . 

Since  y(r,  aQ)  — 1,  we  w that  the  quantity  F(r,on)  introduced  in 
(3.6)  (we  now  make  explicit  its  c endence  on  a ) has  the  property  that 


F(r'V 


1_ 

6 


as  r -*  x 


and  that  (from  the  investigation  f the  ase  L = -1/6  in  Lemma  2)  y(r,  aQ)  — 1 , 
y'(r,  aQ)  — 0 It  :s  Hm  . rhat  the  ext r«  ssion  in  [ . . . ] 

in  (3.  3),  which  we  will  d<  .<?  by  G(t,  u).  has  the  property  that 


as  t — « 
o 


Hence,  given  e > 0,  we  can  ftr.^  r^  sufficiently  large  that 


G(r0’V±  -6  ’ 


and  so,  using  continuous  dependence  on  o,  for  a sufficiently  close  to 
Oq,  we  have 


G(V"'i  ■ 6 * 2£ 


W 


1 * 

L 


From  (3.  3)  we  can  deduce  for  rQ  that 

G(r,  a)  < - ^ + 2e  , (5. 1) 

and  it  is  then  clear  that  y(r,  a)  cannot  become  zero  for  a finite  value  of  r 
or  have  y(°°,  o)  = 0;  for  if  either  of  these  alternatives  were  to  occur,  then 
G(r,  a)  would  become  non-negative  either  at  a finite  value  of  r or  in  the 
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limit  as  r -*•  oo,  and  this  contradicts  (5. 1)  with  £ < 1/12  . The  only  re- 
maining possibility  is  that  y(r,  a)  — 1 as  r — °o  and  the  lemma  is  proved. 

6.  Proof  of  Lemma  5 

We  start  by  making  a change  of  variables  in  (1. 1)  to  take  account  of 
the  fact  that  a is  large.  Thus  we  set 

t = a1/3  r,  Y(t)  = «'2/3  y(r)  , 

and  (1. 1)  becomes 

Y"  + j"  Y1  + {Y  - (a'Z/3  + 1)}Y  = 0 , (6.1) 

where  primes  now  denote  differentiation  with  respect  to  t . The  initial 
conditions  are 

Y(0)  = 0,  Y'  (0)  = 1 , (6.2) 

and  we  consider  in  tandem  with  (6.  I)  the  equation 

YS  f f Yi  * {Yo  - T )Yo  - 0 • <6-  3> 

with  the  same  initial  conditions 

Yq(0)  = 0,  Y^(0)  = 1 . 

If  we  can  show  that  YQ (t)  becomes  negative  for  some  finite  t,  then  since 
(for  sufficiently  large  a)  the  solution  of  (6.  1) -(6.  2)  can  be  made  as  close 
as  we  please  to  YQ(t)  in  any  compact  t-interval,  it  will  follow  that  Y(t) 
also  becomes  negative,  and  the  lemma  is  proved. 
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To  prove  that  becomes  negative,  we  note  that  YQ,  Y^  are 

2 

initially  positive,  and  that  t Y^  is  increasing  until  YQ  meets  the  function 

2 2 2 
2/t  . Thus  Yq  must  meet  2/t  (at  a finite  point),  since  2/t  is  de- 

creasing. 

,2  * , 2 

We  now  ask  whether  Y , after  passing  above  2/t  , meets  2/t 

2 

again.  Suppose  for  contradiction  that  it  does  not.  Then  t Y^  is  decreasing, 

and  it  must  become  negative,  since  otherwise  YQ  is  increasing  and  (6.  3) 

implies  that  Y^  becomes  strictly  negative,  which  ultimately  forces  Y^ 

to  br  negative.  We  thus  have  Y^  positive  and  Y^  negative,  and  so  the 

2 

solution  is  bounded  and  exists  for  all  t,  and  t Y^  -*  L (say)  as  t — oo  , 
with  L < 0 . We  shall  suppose  L finite,  a similar  argument  applying  if 
L = -oo  . Thus  Yq  ~ L/t  , Yq  — M (say),  M finite,  and  it  is  clear  from 
(6.  3)  that  the  only  possible  value  for  M is  M = 0 . Thus 


Y 


0 


~ -L/t 


1 


and  (6.  3)  gives 


(t2Y^)’  „ -L2  , 


Y'  - -L2/e  . 


which  is  a contradiction. 

/ 2 

Thus  Yq  meets  2/t 


a second  time,  and  at  the  point  of  meet  we 


Y-  < -4/t 


3 


must  have 


Now  the  solution  of  the  equation 


2 2 

z"  + — z'  - — z = 0 
1 t2 


with  the  initial  conditions  (at  some  Doint  t ) 

o' 


z('o>  = 2/t*.  z’(tQ)  = -4/t’ 


is  precisely  the  function  2/t  , which  remains  positive  but  tends  to  zero 

2 

as  t — and  comparing  the  equations  for  YQ  and  z = Z/t  , we  have 


(V2>"+f<Vz»'  -T(VZ>  = -^  > 


whence,  as  in  Section  2, 


V2  = MtT  "I  / (‘-4)Yo|u)du 

* *o  1 


with  the  conditions 


<Yo - Z)(V  = °»  • 

This  leads  easily  to  the  inequalities  A < 0,  B > 0,  and  so,  as  t -*  t 


Y - z = At  + — - L(t),  say, 

0 i 

where  L(t)  is  positive  increasing.  Since  A < 0,  this  certainly  implies 


that  Yq  becomes  negative  and  completes  the  proof  of  the  lemma. 
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7 . Proof  of  Lemma  6 

Arguing  as  in  Lemma  4,  we  see  that  the  lemma  will  be  proved  if  we 
can  show  that,  for  a sufficiently  small,  there  is  some  rQ,  possibly  de- 
pending on  a,  for  which 

G(rQ,  a)  < 0 . (7.  1) 

In  order  to  establish  the  existence  of  such  an  r , we  write  (1.  1)  in  the  form 


y"  t-y 


- d+T)y 

r 


2 

-y  , 


(7.2) 


and  note  that  the  solutions  of 


y"  + j y*  - (1  + ~)y  = 0 
r 


er  1 e’r  1 

f(0  = —(1  - -),  g(r)  = — (1  + - 


Using  the  variation  of  constants  formula,  we  can  write  (7.2)  in  the  form 
(taking  account  of  the  initial  conditions) 

y = | or(f  + g)  - \r  f {f(r)  g(t)  - f(t)  g(r)}t2  y2(t)dt  , 


from  which  the  usual  iteration  process  assures  us  that 


3 2 

y - - o(f  + g)  = 0(a  ) as  a - 0 , 


the  estimate  being  uniform  in  any  fixed  interval  [0,  R]  . We  shall  therefore 


certainly  have  established  (7.  1)  if  we  can  show  that  there  exists  a fixed  r 
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(independent  of  a)  for  which  (7.  1)  holds  for  a sufficiently  small  in  the 
form 


G(r0>  a)  < -Kor  , 


where  K is  a positive  constant  independent  of  a and,  in  evaluating  G , 

3 13  3 

we  set  y = -«(f  + g)  . Since  the  term  - y in  G is  0(a  ),  it  will  be 

sufficient  to  prove  that  there  exists  a fixed  rQ  for  which 


2,  * r 

1 ? l ? y (rn)  o 0 , , 

2 y'  (r0)  - 2 7 V -~T~  V / V (t)dt<-Kor 

ro  ° ° 


(7-  3) 


for  a sufficiently  small,  with  y = -o(f  + g).  In  evaluating  the  left-hand 

3 2 

side  of  (7 . 3),  we  can  neglect  the  factor  (—  a)  which  is  common  to  every 

er  l 

term,  and  if  we  replace  f + g by  just  f = — (1  --),  then  we  neglect  only 

terms  which  (for  large  rQ)  are  at  worst  polynomial  in  l/rQ  . With  therefore 
er  1 

y = — (1  - — ),  it  is  standard  to  show  that,  as  r -•>  oo  , 


1 ,2  1 2 y2  2 / 2 3 e2r 

2 7 ~ 2V  ~ 2 + r (t)dt  ~ ‘ 2 4 ’ 


which  certainly  implies  that,  if  r is  fixed  sufficiently  large,  then  the  term 

2r  u 

3 0 4 

arising  from  e /r  on  the  left  of  (7.  3)  will  dominate  any  of  the  other 
terms,  and  so  prove  (7.  3)  and  the  lemma. 
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